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PREFACE 


The  research  work  described  in  this  report  was  conducted  for  the  Defense  Threat  Reduction  Agency 
(DTRA)  under  contract  number  DTRA01-03-D-0014-0030.  The  Contract  Officer’s  Representative  for 
this  effort  was  Mr.  Rick  Fry  of  DTRA’s  Information  Systems  Capability  Development  Office  (J9CBI). 


SUMMARY 


Anatomically  accurate  computational  fluid  dynamics  (CFD)  models  of  the  nasal  passages  of  an  infant  (6 
months  old,  1.3  kg)  and  adult  (7  years  old,  11.9  kg)  rhesus  monkey  were  used  to  predict  deposition  of 
inhaled  nano-  and  micro-particles.  Steady-state,  inspiratory  airflow  simulations  were  conducted  in  each 
model  at  flow  rates  equal  to  100%,  200%,  and  300%  of  the  estimated  minute  volume  for  resting 
breathing.  Particle  transport  and  deposition  simulations  were  conducted  using  the  Lagrangian  method  to 
track  the  motion  of  individual  particles  that  were  passively  released  from  the  nostrils.  Particle  deposition 
fractions  were  higher  in  the  infant  model  than  the  adult  model  at  equivalent  physiological  flow  rates. 
However,  deposition  curves  collapsed  when  accounting  for  differences  in  nasal  geometry  by  plotting 
microparticle  deposition  against  the  Stokes  number  and  nanoparticle  deposition  as  a  function  of  the 
Schmidt  number  and  diffusion  parameter.  Particle  deposition  was  also  quantified  on  major  nasal 
epithelial  types.  Maximum  olfactory  deposition  ranged  from  5-14%  for  1-2  nm  particles  in  the  adult  and 
infant  models,  depending  on  flow  rate.  For  these  particle  sizes,  deposition  on  respiratory/transitional 
epithelia  ranged  from  40-50%.  Increased  deposition  was  also  predicted  for  olfactory  and 
respiratory/transitional  epithelia  for  particle  sizes  >  5  pm  in  the  infant  model  and  >  8  pm  in  the  adult 
model.  Semi-empirical  curves  for  nasal  deposition  were  developed  based  on  the  simulation  data  for 
nano-  and  micro-particle  deposition  to  allow  for  simplified  calculations  of  age -based  nasal  deposition  in 
the  rhesus  monkey  nasal  passages  that  can  be  implemented  into  lung  dosimetry  models. 


1.  INTRODUCTION 


Assessment  of  human  health  risks  from  exposure  to  inhaled  materials  often  relies  on  extrapolation  of 
dose-response  data  from  laboratory  animals.  Nonhuman  primates  are  frequently  used  in  inhalation 
toxicology  studies  due  to  their  similarities  in  respiratory  tract  structure  to  humans  (Desso,  1993;  Hislop 
et  al.,  1984;  Harkema,  1991).  Monkeys  have  been  used  to  assess  effects  from  inhalation  of  reactive 
gases  such  as  ozone,  chlorine  and  formaldehyde  (Carey  et  al.,  2007;  Klonne  et  al.,  1987;  Monticello  et 
al.,  1989),  diesel  exhaust  and  coal  dust  (Nikula  et  al.,  1997),  metal  particles  (Dorman  et  al.,  2006; 
Griffin  et  al.,  1975),  and  infectious  biological  agents  (Day  and  Berendt,  1972;  Druett  et  al.,  1953),  to 
name  a  few  examples.  During  resting  breathing,  the  nasal  passages  serve  as  the  primary  entryway  for 
inhaled  particles  and  gases  and  are  therefore  at  risk  of  insult  from  their  deposition.  An  accurate 
quantification  of  regional  dose  to  the  epithelial  cells  lining  the  nose  is  critical  in  determining  the 
exposure-dose-response  behavior  in  monkeys  and  evaluating  subsequent  risk  for  humans  via 
extrapolation. 

Even  though  monkeys  and  humans  share  many  commonalities,  there  are  differences  in  airway 
morphology  that  affect  the  deposition  of  inhaled  materials.  Many  studies  have  been  conducted  to 
quantify  nasal  particle  deposition  in  humans,  and  several  studies  have  been  conducted  to  determine  how 
differences  in  particle  size  affect  nasal  deposition  and  subsequent  responses  in  monkeys.  Kelly  et  al. 
(2005)  used  a  replica  cast  of  the  nasal  airways  of  an  adult  rhesus  monkey  to  measure  the  deposition 
efficiency  of  particles  in  the  1-10  pm  size  range  and  found  that  deposition  increased  with  flow  rate  and 
particle  size.  Yeh  et  al.  (1992)  also  used  a  replica  cast  to  measure  the  deposition  efficiency  of 
nanoparticles  in  the  5  -  200  nm  size  range.  Nanoparticle  deposition  was  also  measured  in  vivo  by  Yeh  et 
al.  (1997)  in  adult  rhesus  monkeys  and  their  results  agreed  well  with  the  replica  cast  data,  showing  that 
deposition  increased  as  nanoparticle  size  decreased.  Cheng  et  al.  (2008)  measured  respiratory  deposition 
of  2  and  5  pm  droplets  in  cynomolgus  monkeys  and  observed  higher  extrathoracic  airway  deposition  at 
the  larger  particle  size  (58%  versus  39%). 

Particle  size  effects  were  also  noted  following  exposure  to  infectious  agents.  In  the  study  by  Day  and 
Berendt  (1972),  respiratory  tract  responses  in  monkeys  exposed  to  aerosol  particles  containing  tularemia 
(Pasteurella  tularensis)  differed  depending  on  the  particle  size.  Animals  exposed  to  2.1  or  7.5  pm 
particles  displayed  lung  congestion  upon  necropsy,  whereas  effects  in  animals  exposed  to  12.5  or  24  pm 
particles  were  confined  to  the  upper  respiratory  tract.  They  also  found  that  as  particle  size  increased,  an 
increase  in  dose  was  required  to  initiate  a  lethal  infection.  In  the  study  by  Harper  and  Morton  (1953), 
monkeys  exposed  to  aerosolized  spores  of  anthrax  (Bacillus  anthracis)  showed  increased  head  retention 
as  particle  size  increased  from  1  to  12  pm.  Furthermore,  Druett  et  al.  (1953)  reasoned  that  the  site  and 
level  of  infection  was  dependent  on  particle  size  after  exposing  monkeys  to  12  pm  diameter  or  single¬ 
spore  anthrax  particles. 

Experimental  studies  have  provided  invaluable  data  on  the  effects  of  particle  size  on  respiratory 
deposition  characteristics,  yet  each  study  is  limited  by  the  number  of  animals  exposed  and  is  typically 
confined  to  a  narrow  particle  size  range.  Computational  fluid  dynamics  (CFD)  models  allow  for 
simulation  of  airflow  patterns  and  particle  deposition  efficiencies  in  complex  geometries  such  as  those 
found  in  the  upper  respiratory  tract  of  laboratory  animals  and  humans  and  provide  a  valuable 
supplement  to  experimental  work  in  evaluating  dose-response  behavior.  CFD  models  have  been 
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developed  to  study  airflow,  gas  uptake  and  particle  deposition  in  multiple  species,  including  rats,  dogs, 
and  humans  (e.g.,  Ghalati  et  al.,  2012;  Kimbell  et  al.,  1997;  Jiang  and  Zhao,  2010;  Lawson  et  al.,  2012; 
Schroeter  et  al.,  2012).  There  have  also  been  several  CFD  studies  on  vapor  uptake  of  inhaled  chemicals 
in  monkeys.  Kepler  et  al.  (1998)  developed  a  CFD  model  of  the  adult  rhesus  monkey  nasal  passages  to 
study  airflow  and  formaldehyde  uptake.  This  model  was  subsequently  used  by  Kimbell  et  al.  (2001)  to 
compare  formaldehyde  flux  patterns  in  rats,  monkeys,  and  humans.  Corley  et  al.  (2012)  used  a  CFD 
model  of  the  upper  and  lower  respiratory  tract  of  an  infant  rhesus  monkey  to  study  uptake  of  inhaled 
acrolein  vapor.  While  there  have  been  a  couple  efforts  modeling  lung  deposition  of  inhaled  particles  in 
monkeys  (Agharian  et  al.,  2012;  Martonen  et  al.,  2001),  there  have  not  been  any  CFD  studies  to  date  that 
quantify  the  deposition  of  inhaled  particles  in  the  upper  respiratory  tract  of  monkeys.  Nasal  particle 
filtration  is  important  to  understand  to  be  able  to  predict  localized  dose  to  the  nasal  epithelium  and  to 
quantify  the  fraction  of  inhaled  particles  that  reach  the  lung. 

In  this  study,  CFD  models  of  the  nasal  passages  of  an  adult  and  infant  rhesus  monkey  (Macaca  mulatta) 
were  used  to  simulate  inhaled  airflow  and  particle  deposition.  Total  and  epithelial  deposition  fractions 
were  quantified  for  inhaled  nanoparticles  (0.5-1000  nm)  and  microparticles  (1-20  pm).  Semi-empirical 
curves  were  developed  based  on  the  simulation  results  to  extrapolate  nasal  deposition  to  other  exposure 
conditions  and  to  allow  estimation  of  nasal  filtration  for  rhesus  monkey  lung  dosimetry  models. 
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2.  METHODS 


Nasal  airflow  models  were  developed  for  an  adult  (7-year  old  male,  11.9  kg)  and  infant  (6-month  old 
male,  ~1.3  kg)  rhesus  monkey  (Figure  2-1).  The  reconstruction  of  the  adult  model  was  originally 
developed  by  Kepler  et  al.  (1998)  from  serial  Magnetic  Resonance  Imaging  (MRI)  sections  of  the  right 
nasal  passages  (nostril  to  nasopharynx). 


Figure  2-1.  3D  reconstructions  of  the  adult  (top)  and  infant  (bottom)  rhesus  monkey  nasal  passages. 

Approximate  locations  of  olfactory  (dark  grey)  and  squamous  (black)  epithelium  are 
identified.  The  larynx/trachea  region  (dark  grey)  is  identified  in  the  infant  model.  Image  sizes 

are  not  to  scale. 

For  this  study,  the  original  numerical  mesh  developed  by  Kepler  and  colleagues  that  was  comprised  of 
hexahedral  elements  was  removed  and  the  surface  mesh  was  refined  and  smoothed  to  generate  more 
anatomically  realistic  surface  contours  of  the  nasal  passages.  The  infant  rhesus  monkey  model  was 
originally  developed  by  Corley  et  al.  (2012)  and  was  based  on  Computed  Tomography  (CT)  imaging  of 
the  upper  airways  and  MRI  imaging  of  a  lung  cast  to  form  a  contiguous  reconstruction  from  the  nostrils 
through  19  airway  generations  of  the  lung.  For  this  study,  the  model  developed  by  Corley  and  colleagues 
was  truncated  midway  through  the  trachea  to  focus  on  particle  deposition  in  the  upper  respiratory  tract 
(Figure  2-1).  The  major  nasal  epithelial  types  consisting  of  squamous  and  olfactory  epithelia  were 
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mapped  onto  the  surface  mesh  of  each  model  according  to  the  definitions  from  Carey  et  al.  (2007)  and 
Kimbell  et  al.  (2001)  (see  Figure  2-1  and  Table  2-1).  The  remaining  surface  is  composed  of  transitional 
and  respiratory  epithelium. 


Table  2-1.  Geometric  measurements  from  the  monkey  nasal  models 


Region 

Adult  model  Infant  model 

Surface  area  (mm2) 

Whole  nose 

Squamous 

Olfactory 

Larynx/Trachea 

7588.6  2019. 7a 

899.2  204.8 

979.2  331.2 

N/A  458.2 

Volume  (mm3) 

7675.2 

1110.7 

Length  (mm)b 

75.5 

41.4 

Acs  (mm2)c 

142.5 

39.2 

Area  of  nostrils  (mm2) 

88.8 

11.9 

Perimeter  of  nostrils  (mm) 

53.8 

23.8 

Excludes  larynx/trachea  surface  area 

blength  from  the  tip  of  the  nares  to  the  back  of  the  nasopharynx 

ceoronal  cross-sectional  area  of  the  airway  lumen  at  the  beginning  of  the  middle  turbinate 


Numerical  meshes  were  generated  in  each  model  using  ICEM-CFD™  (ANSYS,  Inc.).  Each  mesh 
consisted  of  unstructured  tetrahedral  elements  with  a  thin  3-layer  boundary  of  prism  elements  along 
airway  walls.  The  adult  and  infant  monkey  nasal  models  contained  5.8  and  7.9  million  elements, 
respectively. 

Steady-state,  laminar  inspiratory  airflow  simulations  were  conducted  using  Fluent™  (ANSYS,  Inc.)  to 
solve  the  viscous,  incompressible  Navier-Stokes  equations.  The  segregated  algorithm  was  used  to  solve 
the  continuity  and  momentum  equations;  the  SIMPLEC  method  was  used  for  pressure- velocity 
coupling;  a  second  order  scheme  was  used  to  solve  for  pressure;  and  a  second-order  upwind  scheme  was 
used  to  solve  the  momentum  equations.  Airflow  simulations  were  conducted  with  a  density  of  1.204 
kg/m3  and  a  dynamic  viscosity  of  1.8  x  10-5  kg/m-sec.  Boundary  conditions  in  the  nasal  CFD  models 
consisted  of  a  pre-defined  mass  flow  condition  at  the  nostril(s),  zero  pressure  at  the  outlet,  and  a  no-slip 
condition  on  all  airway  walls.  Airflow  simulations  were  conducted  at  the  minute  volume  (MV),  twice 
the  minute  volume  (2MV),  and  three  times  the  minute  volume  (3MV)  to  span  a  range  of  physiological 
breathing  rates.  Minute  volumes  were  estimated  using  the  allometric  formula  from  Bide  et  al.  (2000): 


MV  =  0.499BW0  809, 
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where  MV  =  minute  volume  (L/min)  and  BW  =  body  weight  (kg).  Using  this  formula,  minute  volumes 
for  the  adult  and  infant  monkeys  were  estimated  to  be  3.7  L/min  and  617  ml/min,  respectively.  Airflow 
simulations  were  therefore  conducted  at  flow  rates  of  1.85,  3.7,  and  5.55  L/min  in  the  adult  model  (flow 
rates  reflect  half  the  whole-nose  flow  rates  since  the  adult  model  consisted  of  only  one  side  of  the  nose) 
and  617,  1234,  and  1851  ml/min  in  the  infant  nasal  model.  Reynolds  numbers  at  the  nostrils  were  307, 
614,  and  921  in  the  adult  model  and  116,  232,  and  348  in  the  infant  model  at  the  respective  flow  rates, 
thereby  justifying  the  use  of  a  laminar  model  for  inhaled  airflow. 

Particle  trajectories  were  calculated  using  the  Discrete  Phase  Model  in  Fluent™  to  solve  the  Lagrangian 
equations  of  motion  which  were  decoupled  from  the  flow  field.  Particle  transport  is  a  function  of  the 
drag  force,  gravity,  and  for  submicron  particles  also  includes  diffusion  effects  due  to  Brownian  motion. 
The  Brownian  force  was  modeled  by  a  Gaussian  white-noise  random  process  model  (Li  and  Ahmadi, 
1992)  and  was  implemented  with  user-defined  functions  in  Fluent  to  improve  upon  the  default  methods 
for  modeling  deposition  of  submicron  particles.  Further  details  of  the  solution  procedures  for  particle 
trajectories  can  be  found  in  Schroeter  et  al.  (2012). 

Uniform  concentrations  of  monodisperse  particles  with  a  density  of  1.0  g/cm3  were  passively  released 
from  the  nostril(s).  Particle  sizes  of  1  nm  to  20  pm  were  simulated.  Nasal  airway  walls  were  defined  to 
be  perfectly  absorbing.  Deposition  fractions  for  the  whole  nose  and  epithelial-specific  (squamous, 
olfactory,  and  respiratory/transitional)  regions  were  computed  by  dividing  the  number  of  particles 
deposited  by  the  number  of  particles  released.  Submicron  particle  deposition  fractions  were  computed 
from  the  average  of  five  independent  simulations. 

Microparticle  deposition  predominantly  occurs  by  inertial  impaction  due  to  the  high  airflow  velocities 
and  narrow  airway  channels  present  in  the  nasal  passages.  In  the  inertial  regime,  deposition  results  can 
be  characterized  by  the  impaction  parameter  (IP),  which  is  a  function  of  aerodynamic  diameter  and  flow 
rate: 


IP  =  pd2Q 


where  p  is  the  particle  density,  d  is  the  particle  diameter,  and  Q  is  the  inhalation  flow  rate.  However,  the 
impaction  parameter  does  not  account  for  different  geometric  conditions  when  comparing  results  from 
models  that  are  similar  in  shape  but  have  different  overall  dimensions.  Instead,  the  Stokes  number  is 
commonly  used  to  compare  inertial  deposition  predictions  from  multiple  models  since  it  includes  terms 
for  the  characteristic  diameter  and  velocity.  The  Stokes  number  can  be  written  as 


Stk  = 


pd2U 
18  judc 


where  U  is  the  characteristic  airflow  velocity,  p  is  the  dynamic  viscosity,  and  d,  is  a  characteristic 
diameter.  Assuming  the  monkey  nasal  passages  retain  their  general  shape  but  grow  in  size  with  age, 
microparticle  deposition  was  compared  in  the  adult  and  infant  monkey  nasal  models  through  the  use  of 
the  Stokes  number  where  the  characteristic  diameter  was  calculated  for  each  model.  Various  forms  of 
the  characteristic  diameter  were  investigated  based  on  geometric  measurements  of  the  models  and  are 
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summarized  in  Table  2-2.  These  included  setting  the  characteristic  diameter  equal  to  the  hydraulic 
diameter  at  the  nostrils,  which  assumes  that  variability  in  deposition  can  be  attributed  to  differences  in 
entrance  effects,  computing  the  diameter  from  the  volume  and  length  scales  of  the  models  (Grgic  et  al., 
2004),  and  defining  the  diameter  in  terms  of  the  surface  area  to  volume  ratio.  We  also  investigated  the 
characteristic  diameter  defined  in  terms  of  the  cross-sectional  area.  Although  Cheng  (2003)  proposed 
using  the  minimal  cross-sectional  area,  a  minimum  value  was  not  apparent  from  our  cross-sectional  area 
calculations  in  the  monkey  nasal  models  (Figure  2-2),  so  the  characteristic  diameter  was  computed  using 
the  cross-sectional  area  at  the  beginning  of  the  middle  turbinate  in  both  models.  All  geometric 
measurements  were  calculated  from  the  numerical  mesh  using  ICEM-CFD™  (ANSYS,  Inc.). 


Table  2-2.  Characteristic  diameters  (mm)  used  in  the  Stokes  number  calculations. 


Characteristic  Diameter 

Adult  model 

Infant  model 

4*An0stril/P  nostril 

6.60 

2.00 

2^VUL 

11.40 

5.85 

4*SAVR1 

4.05 

2.21 

V44s/^ 

13.5 

7.07 

160 

-♦-Adult 


Normalized  distance 

Figure  2-2.  Cross-sectional  areas  in  the  adult  and  infant  rhesus  monkey  nasal  models  as  a  function  of 

normalized  distance  from  the  nostrils  to  the  nasopharynx. 
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For  nanoparticles,  which  primarily  deposit  by  diffusion,  Ingham  (1991)  showed  that  deposition  can  be 
quantified  in  terms  of  the  diffusion  parameter,  A,  and  the  Schmidt  number,  Sc.  The  diffusion  parameter 
is  defined  as: 


Ped 

c  ? 

where  L  is  a  characteristic  length  and  Pe  is  the  Peclet  number,  defined  as 

PeJJL. 

D 

D  is  the  molecular  diffusivity  of  particles  in  air  and  can  be  calculated  from  the  Stokes-Einstein  equation 
(Hinds,  1999).  The  Schmidt  number  is  defined  as 


D 

where  v  is  the  kinematic  viscosity  of  air.  The  characteristic  diameters  defined  above  were  used  to 
calculate  diffusion  parameters  for  comparison  of  nanoparticle  deposition  results  between  the  adult  and 
infant  monkey  models. 
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3.  RESULTS  AND  DISCUSSION 


3.1.  AIRFLOW  SIMULATIONS 

Steady-state  inspiratory  airflow  simulations  were  conducted  in  the  adult  and  infant  monkey  models. 
Airflow  patterns  were  similar  to  the  results  from  Kepler  et  al.  (1998)  and  were  characterized  by  some 
swirling  flow  in  the  vestibule  and  streamlined  flow  in  the  central  nasal  passages  (Figure  3-1).  At  flow 
rates  equal  to  twice  minute  volume,  maximum  airflow  velocities  were  approximately  3-4  m/sec  in  the 
vestibule  and  decreased  to  1-3  m/sec  in  the  central  nose.  The  highest  airflow  velocities  occurred  in  the 
nasopharyngeal  region  where  the  airway  diameters  narrowed  (Figure  3-1).  The  fraction  of  inspiratory 
airflow  directed  to  the  olfactory  region  was  approximately  22%  in  both  models. 


Figure  3-1.  Steady-state  inspiratory  airflow  streamlines  at  twice  minute  volume.  Adult  monkey:  7.4  L/min; 
infant  monkey:  617  ml/min.  Maximum  airflow  velocities  are  5  m/sec. 

The  pressures  at  the  nostril(s)  and  the  nasopharyngeal  outlet  were  computed  to  calculate  the  pressure 
gradient  across  each  nasal  model  under  steady  airflow  conditions  (Figure  3-2).  Pressure  drop  as  a 
function  of  flow  rate  in  the  adult  model  fell  between  the  pressure-flow  curves  from  two  other  studies 
that  measured  pressure  drop  in  adult  monkey  nasal  replica  casts  (Kelly  et  al.,  2005;  Yeh  et  al.,  1992).  At 
equivalent  flow  rates,  pressure  drops  in  the  infant  model  were  significantly  larger  than  the  adult  model 
due  to  narrower  airways  in  the  smaller  model  contributing  more  resistance  to  airflow. 
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Figure  3-2.  Transnasal  pressure  drop  (Pa)  versus  steady  inspiratory  flow  rate  (L/min)  in  the  adult  and 
infant  monkey  nasal  CFD  models  and  the  adult  monkey  nasal  molds  from  Kelly  et  al.  (2005) 

and  Yeh  et  al.  (1990). 

3.2.  NASAL  MICROPARTICLE  DEPOSITION 

Deposition  of  microparticles  in  the  rhesus  monkey  nasal  passages  displayed  sigmoidal  shapes  with 
estimated  deposition  fractions  near  1%  for  1  pm  particles  and  increasing  with  particle  size  (Figure  3-3 
and  Figure  3-4).  At  twice  minute  volume,  nasal  deposition  reached  100%  at  particle  sizes  of  14  pm  and 
9  pm  in  the  adult  and  infant  models,  respectively.  At  equivalent  physiological  flow  rates,  deposition  was 
higher  in  the  infant  model  for  all  particle  sizes.  Higher  nasal  deposition  was  observed  at  higher  flow 
rates  due  to  increased  inertial  effects.  When  deposition  was  plotted  against  the  impaction  parameter, 
deposition  curves  from  multiple  flow  rates  collapsed  for  each  model,  signifying  the  effects  of  inertial 
impaction  on  deposition  (Figure  3-5).  Deposition  values  in  the  adult  model  were  slightly  lower  than  the 
data  from  Kelly  et  al.  (2005).  This  discrepancy  can  partially  be  attributed  to  the  fact  that  the  adult  CFD 
model  used  in  this  study  did  not  contain  a  larynx  or  trachea  which  would  result  in  less  deposition,  and 
may  also  be  due  to  inter-subject  variability  between  the  two  models.  Deposition  in  the  infant  model  was 
significantly  higher  than  the  data  from  Kelly  et  al.  (2005)  due  to  the  smaller  model  size. 


2  4  6  8  10 

Flow  rate  (L/min) 


9 


0  5  10  15  2 

Particle  size  (\irr\) 


Figure  3-3.  Microparticle  nasal  deposition  predictions  in  the  monkey  nasal  models  as  a  function  of  particle 

size  -  adult  model. 


Figure  3-4.  Microparticle  nasal  deposition  predictions  in  the  monkey  nasal  models  as  a  function  of  particle 

size  -  infant  model. 
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Figure  3-5.  Microparticle  nasal  deposition  predictions  in  the  adult  and  infant  monkey  nasal  models  as  a 
function  of  the  impaction  parameter.  Deposition  data  from  an  adult  rhesus  monkey  nasal  cast 

is  included  (Kelly  et  al.,  2005). 


Nasal  deposition  fractions  in  the  adult  and  infant  models  were  also  plotted  against  the  Stokes  number 
using  various  forms  of  the  characteristic  diameter  from  Table  2-2  (Figure  3-6,  Figure  3-7,  and  Figure 
3-8,  and  Figure  3-9).  Defining  the  characteristic  diameter  to  be  equal  to  the  hydraulic  diameter  at  the 
nostrils  did  not  fully  account  for  differences  in  deposition  between  the  two  models  (Figure  3-6). 
Although  nostril  inlet  conditions  certainly  influence  deposition,  this  form  of  the  characteristic  diameter 
did  not  account  for  any  downstream  differences  in  model  geometry  that  could  also  affect  the  deposition 
of  particles  throughout  the  nose.  Use  of  other  forms  of  the  characteristic  diameter  resulted  in  a  better 
collapse  of  the  deposition  curves  when  plotted  against  the  Stokes  number  and  were  able  to  account  for 
differences  in  model  size  on  deposition  (Figure  3-6,  Figure  3-7,  and  Figure  3-8).  Similar  findings  were 
observed  by  Grgic  et  al.  (2004)  and  Storey-Bishoff  et  al.  (2008)  when  accounting  for  inter-subject 
variability  in  human  oral  airway  deposition  and  infant  nasal  deposition,  respectively.  Storey-Bishoff  et 
al.  (2008)  compared  numerous  forms  of  the  characteristic  diameter  and  found  that  the  volume-surface 
ratio  definition  provided  the  best  collapse  of  infant  nasal  deposition  data.  The  characteristic  diameter  dc 
=  4*SAVR_1  (SAVR  =  surface  area  to  volume  ratio)  was  used  in  further  calculations  since  a  good 
collapse  of  the  data  was  observed  and  since  nasal  volume  and  surface  area  can  be  estimated  from  body 
weight  and  can  therefore  be  used  to  extrapolate  results  to  other  monkey  ages  (Asgharian  et  al.,  2012). 
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Figure  3-6.  Total  nasal  deposition  results  of  inhaled  microparticles  in  the  adult  and  infant  monkey  nasal 
models  as  a  function  of  Stokes  number.  Various  forms  of  the  hydraulic  diameter  were  used  in 

Stokes  number  calculations:  4*An0strii/Pnostrii 
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Figure  3-7.  Total  nasal  deposition  results  of  inhaled  microparticles  in  the  adult  and  infant  monkey  nasal 
models  as  a  function  of  Stokes  number.  Various  forms  of  the  hydraulic  diameter  were  used  in 

Stokes  number  calculations:  2^/v  /  tzL 
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Figure  3-8.  Total  nasal  deposition  results  of  inhaled  microparticles  in  the  adult  and  infant  monkey  nasal 
models  as  a  function  of  Stokes  number.  Various  forms  of  the  hydraulic  diameter  were  used  in 

Stokes  number  calculations:  4*SAVR  1 
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Figure  3-9.  Total  nasal  deposition  results  of  inhaled  microparticles  in  the  adult  and  infant  monkey  nasal 
models  as  a  function  of  Stokes  number.  Various  forms  of  the  hydraulic  diameter  were  used  in 

Stokes  number  calculations:  /  n 
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3.3.  NASAL  NANOPARTICLE  DEPOSITION 


Nanoparticle  simulation  results  were  compared  with  the  experimental  data  of  Yeh  et  al.  (1992),  where 
deposition  fractions  of  1.7  -  100  nm  particles  were  measured  in  a  replica  cast  of  the  nasal  airways  of  an 
adult  rhesus  monkey  (body  weight  8.5  kg).  At  steady  inhalation  flow  rates  of  1.75,  3.5,  and  7.0  L/min, 
Yeh  et  al.  (1992)  measured  nasal  deposition  of  1.7  nm  particles  to  range  from  78-87%  (Figure  3-10). 
Deposition  decreased  as  particle  size  increased  due  to  the  diminished  effects  of  diffusion.  Deposition  of 
100  nm  particles  was  <  11%  at  all  flow  rates.  CFD  simulations  of  nanoparticle  deposition  in  the  adult 
rhesus  monkey  model  at  flow  rates  of  1.85,  3.7,  and  7.4  L/min  compared  well  with  the  experimental  data 
(Figure  3-10).  A  slight  under-prediction  was  observed  for  particle  sizes  <  10  nm.  This  can  be  partially 
attributed  to  the  slightly  higher  flow  rates  used  in  the  CFD  model  and  the  fact  that  the  nasal  replica  cast 
used  by  Yeh  et  al.  (1992)  also  included  the  pharynx  and  larynx. 

^  1  ♦  Yeh  et  al.  (1992):  7.0  L/min 


0.9  \  o  Yeh  et  al.  (1992):  3.5  L/min 

ss 
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Figure  3-10.  Comparison  of  nanoparticle  deposition  predictions  in  the  adult  rhesus  monkey  nasal  CFD 

model  with  the  experimental  measurements  of  Yeh  et  al.  (1992). 

CFD  simulations  of  nanoparticle  deposition  were  conducted  in  the  adult  and  infant  rhesus  monkey 
models  at  flow  rates  equal  to  the  MV,  2MV,  and  3MV  for  particle  sizes  between  0.0005  and  1  pm  (0.5  - 
1000  nm)  (Figure  3-11).  At  equivalent  physiological  flow  rates,  nanoparticle  deposition  in  the  infant 
model  was  greater  than  the  adult  model.  For  each  model,  deposition  fractions  were  lower  at  higher  flow 
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rates  due  to  decreased  residence  times  of  nanoparticles  in  the  nasal  passages.  Deposition  of  0.5  nm 
particles  was  predicted  to  be  100%  in  the  infant  model  at  all  flow  rates.  In  the  adult  model,  0.5  nm 
deposition  ranged  from  89-99%  depending  on  flow  rate.  Flow  rate  differences  were  most  evident  for 
particle  sizes  <10  nm.  Deposition  of  1  nm  particles  ranged  from  65-87%  in  the  adult  model  and  from 
91-99%  in  the  infant  model.  Deposition  decreased  dramatically  as  particle  size  increased  to  from  1  nm 
to  10  nm.  Deposition  of  10  nm  particles  ranged  from  7-16%  in  the  adult  model  and  from  15-28%  in  the 
infant  model.  As  particle  size  increased  above  10  nm,  deposition  slowly  decreased  further  and  was  <  4% 
for  particle  sizes  >  100  nm  at  all  flow  rates. 


0.0001  0.001  0.01  0.1  1 

Particle  Size  (nm) 

Figure  3-11.  Nanoparticle  nasal  deposition  predictions  in  the  adult  and  infant  monkey  nasal  CFD  models. 

To  further  compare  nanoparticle  deposition  in  the  two  models,  deposition  fractions  were  plotted  against 

the  term  Sc“ Afc,  where  a  and  b  are  fitted  constants.  Variability  in  nanoparticle  deposition  was  reduced 
considerably  with  the  constants  a  =  -0.08,  b  =  0.25  (Figure  3-12).  In  this  example,  the  characteristic 
diameter  in  the  diffusion  parameter  calculation  was  defined  as  dc  =  4*SAVR  1  to  be  consistent  with  the 
definition  in  the  microparticle  deposition  results.  Because  airway  dimensions  were  included  in  the 
definitions  of  the  Schmidt  number  and  diffusion  parameter,  this  relationship  should  apply  to  a 
population,  assuming  that  the  shape  of  the  nasal  passages  remains  the  same  but  gets  larger  with  age. 
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Figure  3-12.  Nanoparticle  nasal  deposition  predictions  in  the  adult  and  infant  monkey  nasal  CFD  models 
as  a  function  of  the  diffusion  parameter,  d.  The  characteristic  diameter  was  dc  =  2tJv/ttL. 

3.4.  REGIONAL  PARTICLE  DEPOSITION 

Micro-  and  nano-particle  deposition  fractions  were  tabulated  for  the  major  nasal  epithelial  types: 
squamous,  olfactory,  and  respiratory/transitional  epithelia.  Deposition  was  also  quantified  in  the 
larynx/trachea  region  of  the  infant  model.  Deposition  curves  for  the  entire  particle  size  range  considered 
in  this  study  (0.5  nm  to  20  pm)  are  shown  in  Figures  3-13  to  3-21.  Total  nasal  deposition  curves 
displayed  a  classic  U-shape  that  is  seen  in  other  species,  with  high  deposition  for  small  nanoparticles 
and  large  microparticles  (Figure  3-13  and  Figure  3-14).  Deposition  was  very  low  (<  4%)  for  the  fine 
particle  size  range  of  100  nm  to  1  pm.  As  particle  size  increased  above  1  pm,  deposition  fractions 
increased  sharply.  Similar  behavior  was  observed  on  squamous  epithelia,  with  moderately  high 
deposition  for  small  nanoparticles  and  high  deposition  for  large  microparticles  (Figure  3-15  and  Figure 
3-16).  Deposition  of  1  nm  particles  was  predicted  to  range  from  16-36%  and  from  24-43%  in  the 
squamous  region  in  the  adult  and  infant  models  respectively,  depending  on  flow  rate.  Nanoparticles 
diffuse  to  airway  walls  upon  entering  the  nasal  passages,  which  leads  to  enhanced  deposition  in  the 
anterior  section  of  the  nose  that  is  covered  by  squamous  epithelia.  Large  particles  (>  5  pm)  experience 
deposition  by  inertial  impaction  in  this  region  due  to  high  airflow  velocities  and  rapidly  changing  nasal 
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shapes.  Squamous  deposition  was  higher  in  the  infant  model  than  in  the  adult  model  at  equivalent 
physiological  flow  rates. 


0.0001  0.001  0.01  0.1  1  10 
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Figure  3-13.  Total  nasal  deposition  predictions  in  the  monkey  nasal  CFD  models:  adult. 


Figure  3-14.  Total  nasal  deposition  predictions  in  the  monkey  nasal  CFD  models:  infant. 
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Figure  3-15.  Deposition  predictions  on  squamous  epithelium  in  the  monkey  nasal  passages:  adult. 
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Figure  3-16.  Deposition  predictions  on  squamous  epithelium  in  the  monkey  nasal  passages:  infant. 
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Deposition  on  olfactory  epithelium  is  important  to  quantify  since  there  is  concern  for  translocation  of 
nanoparticles  and  certain  metals  along  olfactory  neurons  to  the  olfactory  bulb  (Dorman  et  al.,  2006; 
Lewis  et  al.,  1994;  Oberdoster  et  al.,  2004).  In  addition,  drugs  depositing  in  the  olfactory  region  can  be 
transported  to  the  central  nervous  system  through  the  olfactory  pathway  (Thorne  et  al.,  2008).  For  the 
particle  sizes  and  flow  rates  considered  in  this  study,  deposition  on  olfactory  epithelium  was  very  low  (< 
14%).  Olfactory  deposition  displayed  a  bimodal  behavior  with  higher  deposition  occurring  at  very  small 
nanoparticle  sizes  (<  3  nm)  and  very  large  microparticle  sizes  (>  10  pm)  (Figure  3-17  and  Figure  3-18). 
The  highest  olfactory  deposition  in  the  adult  model  of  8%  occurred  at  a  particle  size  of  1  nm  at  the 
highest  flow  rate.  Nanoparticle  olfactory  deposition  was  slightly  higher  in  the  infant  model  with  a  peak 
of  14%  occurring  at  a  particle  size  of  1  nm  at  the  highest  flow  rate.  Similar  behavior  was  predicted  for 
particle  deposition  in  the  rat  nasal  passages  with  the  highest  olfactory  deposition  occurring  in  the  1-3  nm 
size  range  (Garcia  and  Kimbell,  2009;  Schroeter  et  al.,  2012).  However,  this  does  contrast  with  predicted 
olfactory  nanoparticle  deposition  in  humans,  which  was  estimated  to  be  <  1%  (Shi  et  al.,  2006). 
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Figure  3-17.  Deposition  predictions  on  olfactory  epithelium  in  the  monkey  nasal  passages:  adult. 
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Figure  3-18.  Deposition  predictions  on  olfactory  epithelium  in  the  monkey  nasal  passages:  infant. 

Deposition  on  the  remaining  nasal  surface  comprised  of  respiratory  and  transitional  epithelia  also 
displayed  a  bimodal  behavior,  with  peak  deposition  occurring  at  particle  sizes  of  0.5-2  nm  and  in  the  7- 
20  pm  range  (Figure  3-19  and  Figure  3-20).  Deposition  in  this  region  was  much  higher  than  in  the 
olfactory  region  due  to  the  large  surface  area  and  higher  proportion  of  airflow  passing  over  these 
epithelial  types. 
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Figure  3-19.  Deposition  predictions  on  respiratory/transitional  epithelium  in  the  monkey  nasal  passages: 

adult. 


Figure  3-20.  Deposition  predictions  on  respiratory/transitional  epithelium  in  the  monkey  nasal  passages: 

infant. 

Nanoparticle  deposition  in  the  larynx/trachea  region  of  the  infant  model  was  <  12%  and  reached  a 
maximum  at  a  particle  size  of  2  nm  (Figure  3-21).  Peak  microparticle  deposition  was  significantly 
higher  at  35%,  occurring  at  the  highest  flow  rate  at  a  particle  size  of  7  pm.  Microparticle  deposition  was 
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elevated  in  the  larynx  area  due  to  the  airway  constriction  resulting  in  higher  airflow  velocities  in  this 
region,  which  increases  deposition  due  to  impaction. 
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Figure  3-21.  Deposition  predictions  in  the  larynx/trachea  region  in  the  infant  monkey  nasal  CFD  model. 
3.5.  SEMI-EMPIRICAL  NASAL  DEPOSITION  CURVES 

The  CFD  simulation  data  was  used  to  develop  semi-empirical  curves  for  total  nasal  deposition.  Particle 
deposition  in  the  upper  airways  of  the  respiratory  tract  directly  affects  the  amount  that  is  available  for 
deposition  in  the  lung.  It  is  desirable  to  construct  analytical  relationships  for  the  predicted  deposition 
fraction  in  the  nasal  passages  as  functions  of  airflow  and  particle  variables  that  can  be  implemented  in 
lung  deposition  models  for  rhesus  monkeys  (Asgharian  et  al.,  2012).  Two  such  relationships  are  needed: 
one  for  sub-micrometer-sized  particles  for  which  deposition  occurs  mainly  by  Brownian  diffusion,  and 
also  for  fine  and  coarse  particles  which  deposit  by  inertial  impaction  in  the  nasal  passages.  Deposition 
by  sedimentation  is  negligible  in  the  upper  airways  because  of  small  particle  residence  times. 

For  micrometer-sized  particles,  nasal  deposition  curves  collapsed  when  plotted  versus  the  Stokes 
number  with  the  characteristic  diameter  defined  in  terms  of  the  surface  area-volume  ratio  (Figure  3-7). 
Likewise,  nanoparticle  deposition  collapsed  when  plotted  as  a  function  of  the  Schmidt  number  and 
diffusion  parameter  (Figure  3-12).  For  consistency  with  the  lung  geometry  model  and  ease  of 
implementation,  deposition  efficiencies  can  also  be  plotted  versus  the  impaction  parameter,  pd2Q,  and  in 
terms  of  the  diffusion  coefficient  (D)  and  flow  rate  (Q)  for  microparticles  and  nanoparticles, 
respectively.  Predicted  deposition  efficiencies  of  inhaled  particles  in  the  infant  nasal  passages  (which 
includes  deposition  in  the  larynx)  were  fit  to  sigmoidal  functions  as  given  below: 
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— Infant:  2MV 

—  Infant:  3MV 
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where  r|i  and  rjd  represent  deposition  efficiencies  by  inertial  impaction  and  diffusion,  respectively 
(Figure  3-22  and  Figure  3-23).  The  impaction  parameter,  pd2Q,  in  the  first  equation  is  given  in  units  of 
g-pm2/sec  while  the  flow  rate,  Q,  and  diffusion  coefficient,  D,  in  the  second  equation  are  given  in 
cmVsec  and  cm2/sec,  respectively.  In  addition,  optimized  parameter  values  for  the  coefficients  are 
shown  in  Table  3-1. 


pd2Q,  cm3/s 


Figure  3-22.  Empirical  curve  for  nasal  deposition  of  microparticles  as  a  function  of  the  impaction 

parameter.  Data  points  represent  CFD  simulations  results  from  the  infant  monkey  model. 
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Figure  3-23.  Empirical  curves  for  nasal  deposition  of  nanoparticles  as  a  function  of  the  diffusion 

parameter.  Data  points  represent  CFD  simulation  results  from  the  infant  monkey  model. 


Table  3-1.  Optimized  parameter  values  for  the  sigmoidal  functions  for  nasal  deposition  by  inertial 

impaction  (ip)  and  diffusion  (pd). 


Parameter 

Value 

Impaction  (ip) 

Ai 

0.863 

Bi 

441.540 

Ci 

3.289 

Diffusion  (r|d) 

Fi 

0.969 

Gi 

2549.821 

Hi 

0.901 

Ki 

-0.958 

Assuming  that  the  monkey  nasal  passages  retain  their  shape  during  growth  and  only  change  in  size,  the 
above  deposition  efficiency  formulas  may  be  extended  to  other  ages  by  correcting  for  the  airway  volume 
and  surface  areas  (Asgharian  et  al.,  2012)  to  reflect  the  following  change  in  parameter  values: 
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where  Si,  Vi  and  S2,  V2  correspond  to  the  nasal  surface  area  and  volume  of  the  infant  monkey  model 
(see  Table  2-1)  and  a  monkey  of  different  age,  respectively. 
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4.  CONCLUSIONS 


CFD  models  of  the  nasal  passages  of  an  infant  and  adult  rhesus  monkey  were  used  to  predict 
deposition  of  inhaled  nanoparticles  and  microparticles.  Nasal  deposition  was  very  different 
between  the  two  models  due  to  their  size  differences.  However,  when  deposition  fractions  were 
plotted  against  appropriate  non-dimensional  parameters,  the  variability  between  models  was 
substantially  reduced.  Although  only  two  models  were  used  in  this  analysis,  these  results  suggest 
that  nasal  deposition  curves  can  be  extrapolated  to  rhesus  monkeys  of  other  ages  assuming  that 
the  nasal  passages  retain  their  shape  with  age  and  that  their  nasal  surface  areas  and  volumes  can 
be  measured  or  estimated.  There  may  be  significant  inter-individual  variability  between  monkeys 
that  can  be  explored  once  other  nasal  geometries  become  available. 

The  semi-empirical  nasal  filtration  equations  were  implemented  into  the  Multiple-Path  Particle 
Dosimetry  (MPPD)  lung  model  that  was  recently  developed  for  rhesus  monkeys  (Asgharian  et 
al.,  2012).  This  model  can  be  downloaded  at  the  following  address: 

http://www.ara.com/products/mppd.htm 
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DEFINITIONS,  ACRONYMS,  AND  ABBREVIATIONS 


BW 

Body  Weight 

CFD 

Computational  Fluid  Dynamics 

CT 

Computed  Tomography 

DTRA 

Defense  Threat  Reduction  Agency 

Fluent™ 

CFD  software  sold  by  ANSYS,  Inc. 

ICEM-CFD™ 

Mesh  generation  software  sold  by  ANSYS,  Inc. 

IP 

Impaction  Parameter 

MPPD 

Multiple  Path  Particle  Dosimetry 

MRI 

Magnetic  Resonance  Imaging 

MV 

Minute  Volume 

SAVR 

Surface  Area  to  Volume  Ratio 

SIMPLEC 

Method  in  the  Fluent™  software 

2MV 

Twice  minute  volume 

3MV 

Three  time  the  minute  volume 
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